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ABSTRACT 

The space mission PLATO will usher in a new era of exoplanetary science by ex¬ 
panding our current inventory of transiting systems and constraining host star ages, 
which are currently highly uncertain. This capability might allow PLATO to detect 
changes in planetary system architecture with time, particularly because planetary 
scattering due to Lagrange instability may be triggered long after the system was 
formed. Here, we utilize previously published instability timescale prescriptions to de¬ 
termine PLATO’s capability to detect a trend of decreasing planet frequency with age 
for systems with equal-mass planets. For two-planet systems, our results demonstrate 
that PLATO may detect a trend for planet masses which are at least as massive as 
super-Earths. For systems with three or more planets, we link their initial compactness 
to potentially detectable frequency trends in order to aid future investigations when 
these populations will be better characterized. 
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1 INTRODUCTION 

Nearly all exoplanets orbit stars whose ages are poorly con¬ 
strained. This situation is unfortunate because the potential 
to learn about the dynamical evolution of planetary systems 
through the host star’s evolution has yet to be realized. Ac¬ 
curate stellar ages may provide key constraints and insights 
on the formation and fate of planets. 

Tidally-influenced, hot exoplanets highlight the impor¬ 
tance of this link. Hundreds of confirmed planets have orbital 
periods less than 10 days, where tidal interactions between 
the planet and star can affect the planet’s o rbital parame - 
ters in a myriad of ways (e.g. see Fig. 4 of IOgilviell2014lf . 
Because the timescales for shrinking an orbit are sensitively 
dependent on the orbital period, accurate stellar ages would 
crucially constrain the dynamical histories of these systems. 
Further, the detai ls of tidal theory are controv ersial and re¬ 
cent analyses ('e.g. lEfroimskv fe Makarovll2013l ') demonstrate 
the dangers of adopting the heretofore widely-used constant 
quality factor. Adding stellar age to the list of known system 
parameters would enable greater resolution on this topic. 

However, the insights obtained from accurate stel¬ 
lar ages are not restricted to substellar companions with 
such close orbits. Single planets or brown dwarfs which 
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are out of the reach of stellar tides (with semima¬ 
jor axes beyond about 0.1 au) but within the ellip¬ 
soid of the gravitational influence of the parent star 
within the Milky Way (within about 10 5 au) could at¬ 
tain currently-observed orbital properties from (i) a relic 
planet formation process such as core a ccretion, gravi¬ 
tational instability or d isc migration (see iBaruteau et al.l 
120131 : [Helled et aI1l2013l . for recent reviews), (ii) a relic star 
formatio n process such as star-plan e t scattering in birth 
clusters JPfahl fc Muters paughl [2 0061: IStmrzem et ahl 120091: 


Malmberg et alj ~ 201ll : Parker fc Quana 20121 : Hao et al.l 

2013T) . (iii) gravitational interaction with other planets 


durin g a star’s main sequence lifetime (see iDavies et al.1 
l20l4 for a recent review) or /and (iv) perturbation s from 
Galactic field star flybys fZ akamska & Tremaine 200J; 
Fregeau et al.l 120061 : iBolev et al.l 1201 21: lYeras fc Moeckel 


2012) and Galactic tides llHeisler fc Tremainel 19861 : Brassei 


2001 ; IVeras fc EvansH20lJ T 


All these processes are time-dependent and hence could 
benefit from accurate stellar ages. Techniques are being de - 
veloped to constrain stellar ages (ILebreton fc GoupilH2014r ) 
to within 10 per cent of a typical main sequence lifetime 
of a Solar-mass star by utilizing astroseismology. Although 
these techniques can be applied to high-precision p hotomet- 
ric d ata taken from the C 0 R 0 T jBag lin et al1 l2002l i and Ke¬ 
pler dBorucki et al.ll2010l : iKoch et al.l "oid f space missions, 
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the number of suitable stars from these missions is just a 
handful. 

The PLATO space mission (iR.auer et al.ll201 4i). due for 
launch in 2024, will achieve similarly tight age constraints, 
but for a much larger sample of stars. One of the main aims 
of PLATO is to achieve 10% precision in the age of 20,000 
main sequence stars of spectral types F5-K7. In order to 
achieve this goal, PLATO will produce internal stellar mass 
distributions for different stars by combining radius mea¬ 
surements derived from Gaia data with the oscillation fre¬ 
quencies observed in the PLATO light curves. 

Other techniques which will be used to derive stellar 
ages, such as gyrochronology, will be assessed and calibrated 
from this “primary” sample and then applied to the hun¬ 
dreds of thousands of main sequence stars that PLATO will 
observe during the course of its 6 year mission (albeit at a 
lower age precision). Moreover, PLATO will, for the primary 
sample, allow for direct testing of stellar evolution tracks. 

Such a large sample will enable us to detect trends in 
exoplanetary science that remain currently hidden in the 
noise. One of these trends is frequency of planets versus 
main sequence age. Not all systems with multiple planets 
will remain stable over their main sequence lifetimes. In our 
own Solar system, Mercury has a one to few per cent chance 
of dest abiliz ing the inner planets (lLaskar fe Gastineaull20oH ; 
IZeebel l2015h. Other sy s tems are closer to instability; Sec¬ 
tion 6.3 of I Veras et ahl (l2013h discuss specific examples for 
planets discovered by Doppler radial velocity with semima¬ 
jor axes that exceed 1 au. Overall, the frequency of planets 
might decrease with time. 

The trend for planets detected by transit photometry 
within 1 au but beyond the limit at which tidal forces 
cease to become significant (at about 0.1 au) should not 
be qualitatively different. On 9 Jul 2015, the Exoplanets 
Data Explorer (http://exoplanets.org) listed 509 substel- 
lar objects in this semimajor axis range t hat are compo¬ 
nents of m ultiple-planet systems. Ov erall, iFang fe Margotl 
(l2012fl and I Ballard fe Johnsonl (120141 ') both estimate that 
about half of Kepler systems are multiple-planet systems, 
the latter considering only M dwarfs. Seven examples of po¬ 
tential planets in multi-planet systems are the seven can¬ 
didate planets of KIC 11442793, all of which transit be- 
tween 0.074 au and 1 .0 au but do not have measured masses 
(ICabrera et alj(20141 '). In fact, those authors conclude that 
the two outermost candidates are indeed planets due to dy¬ 
namical stability constraints on their masses, highlighting 
how close to instability some of these systems might reside. 

In this brief paper, we combine stability prescriptions 
from previous investigations with the detection capabilities 
of the PLATO mission. Section 2 identifies these prescrip¬ 
tions and Section 3 performs the analysis. We conclude in 
Section 4. 


2 STABILITY OF PLANETARY SYSTEMS 

Two-planet systems represent a robust starting point with 
which to study decreasing planet frequency with time along 
the main sequence. Observers have identified hundreds of 
these systems, and theorists have established analytical con¬ 
straints which would otherwise be unattainable with addi¬ 
tional planets. 


2.1 Two-planet stability 


A two-planet system is a generally-unso lv able three- 
body problem. Howev er, iMarchal fe Saaril (119751 1 and 
iMarchal fe Boziil d 1982fl identified allowable regions of mo¬ 
tion in this problem based on the angular momentum and 
energy of the system. This categorization proved partic¬ 
ularly v aluable after the confide nt discovery of extrasolar 
planets dWolszczan fe FraiJll992|j for assessi ng their stabil¬ 
ity. M any subsequent studies, summarized in I Geor gakarakoa 
(1200811 . considered relevant special cases and extensions. 

The allowable regions of motions are partially parti¬ 
tioned by the concept of Hill stability. Two planets are said 
to be Hill stable if their orbits never cross. This requirement 
does not preclude the inner planet from crashing into the 
star, nor the outer planet from escaping from the system. If 
the planets do remain ordered and do not suffer a collision 
nor ejection, then the system is said to be Lagrange stable. In 
Hill unstable cases, when the initial orbits do cross, then the 
outcome is uncertain: the system could destabilize or alter¬ 
natively the planets might remain inside of a mean motion 
resonance (like Neptune and Pluto). This latter possibility 
presupposes the planets were originally captured into a res¬ 
onance, most likely due to protoplanetarv disc dissipation, 
but not necessarily (iRavmond et al.ll2008l l. 

Nevertheless, in most cases (as a fraction of available 
parameter space), planets which cross orbits will become vi¬ 
olently unstable. In this respect, the critical separation 
of the initial planetary orbits that determines the Hill sta¬ 
bility boundary has occasionally been t reated, erroneo usly, 
as a type of global stability boundary. iMarzar 1 (l2014l! has 
demonstrated why Hill stability should not be treated as a 
global boundary through a frequency map analysis. In some 
cases though, the Hill stability boundary may be considered 
a dividing line between “quick instability” and “slow insta¬ 
bility”. This notion may be carried further: one can then de¬ 
fine a Lagrange stability boundary Al > Ah beyond which 
Lagrange instability does not occur, or occurs on a timescale 
which exceeds the main sequence stellar lifetime. 

The lack of quantification in these statements re¬ 
sults from a dearth of dedicated studies of Lagrange in¬ 
stability. We summ arize previous i nvestigati o ns as fol¬ 
lows: Chapter 11 of iMarchall (Il990l j. lAnosove] (Il996l l and 
Li et all d2009l) focused on t he esc ape of the outer body; 


Kholshevnikov feJKuznctsov teonT) considered the mass 
dependence: iDeck et alj ~i 2013li ~ linked Lagrange instabil¬ 


ity (albeit defined slightly different ly ) with mean mo¬ 
tion r esonance overlap; an d Barnes_&_Greenberg (2006, 
2007li . IRavmond et al] (l2009l'l. Koppara pu fe Barnes (I201C 


Veras et al.l ( 2013ll and Veras fe Musti ll (2013) evaluated 


the extent of Lag range inst ability beyond the Hill stabil¬ 
ity limit. Recently IPetrovichl (1201511 unveiled a fairly general 
Lagrange instability criterion for two unequal-mass, highly 
eccentric and nearly coplanar planets that performs signifi¬ 
cantly bett er than previous widely-used criteria from stellar 
dynamics (lEggleton fe Kiseleval Il995l : iMardling fe Aarsethl 
l200lll. Although exc ellent analytic approximations to Ah 
exist (iDonniso nl201lll . the same is not true for Al. One must 


1 The boundary of the Hil l radius of a single pla net is a related 
concept (see Appendix B of Pearce &; Wvattj|2014r) , but not equiv¬ 
alent to Ah- 
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find Al empirically through simulations. However, simula¬ 
tions have revealed that this Lagrange stability boundary is 
neb ulous. _ 

IVeras fc Mustilll (l2013i ) found that although Al is diffi¬ 
cult to constrain, the minimum time to Lagrange instability 
is well-established empirically as a function of planet mass 
only, for equal-mass planets. They defined the Lagrange in¬ 
stability timescale as the time at which either the inner 
planet collided with the stellar surface or when the outer 
planet’s orbit became hyperbolic. If x U ns represents the min¬ 
imum number of initial orbits of the inner planet before the 
onset of Lagrange instability, then 


/ \ — 0 . 181 ± 0.003 

log 10 (z una ) ~ (5.23 ± 0.04) m // Mq J > W 

where /x is the planet/star mass ratio, for both planets, and 
Mj is the mass of Jupiter. This formula appears to hold 
true for any separation in-between Ah and Al, for, at least, 
orbits with initial eccentricities less than 0.3, and is valid for 
both planets and brown dwarfs. The equation is based on 
a fit down to masses of O.lMj, and will be used directly in 
conjunction with PLATO sensitivities in Section 3. 


2.2 Stability of three or more planets 


In systems with three or more planets, there is no known 
analytical constraint on Hill stability, and the distinction 
between Hill instability and Lagrange instability is rarely 
asserted. Instead, researchers have relied on empirical fits 
to A-body simulations to help characterize the potential fu¬ 
ture stability of a system. Analytical motivation for estab- 
lishi ng these fits is also rare. For an exception, see Section 
3 of IZhou et al] (120071 1. who demonstrate how an analyti¬ 
cal treatment from first principles can reproduce the same 
scalings obtained from their empirical fits. 

The functional form predominately chosen for the global 
instability fit for thre e or more planets is equivalent to 
dChambers et ah 19961 : ISmith fe Lissaueij 120091 ; iFunk et all 
120101 : IPu fc Wull2015l f 


logic (* uns ) ~ bp + c, 


( 2 ) 


Zhou et al.l 

20071: IChatteriee et al. 

Mustill et al 

201411. In equation ( 


Faber&Quillc 


iUenj 

illenl 


2007 


2011 


ted constants and /? refers to the number of mutual Hill 
radii by which consecutive pairs of planets are sepa¬ 
rated. A mutual Hill radius is a distance that is not de- 
fined uniformly throu gho ut the literature (contrast ; e.g., 
IChambers et al.l 19961 and Marzari fc Weidenschillina 120021 
with Smith fe Lissaue H |2009l, and see the discussion at the 
end of Section 2 of Mustill et al . IH). Regardless, ef¬ 
fectively, /3 determines the extent of packing in a plane¬ 
tary system dBarnes fe Ravmondll2004l; iBarnes et all 120081 : 


iRavmond et all 20091 : Veras fc Ghnsicka 20151 ). and can be 

expressed through masses and semimajor axes. Although the 
exact definiti on of / j varies (s ee also the sentence containing 
equation 5 in IPu fe Wull2015l b overall if /3 is too small, then 
the system will become unstable almost immediately. Alter¬ 
natively, if /3 is too large, the system will not experience 
instability during the star’s main sequence lifetime. 


Like Equation ©, Equation ([2]) is applicable only for 
equal-mass planets. However, Equation ([2]) is more versatile, 
having been applied to systems with for example 3, 5, 10 
and even 20 planets. The unequal-planet mass case greatly 
increases the available phase space to explore, and has so- 
far largely been ignored (we do so here as well). Although 
Equation m may be more versatile, this general form does 
not appear to constrain stability as well as Equation m, 
because the latter is for the specific case of two equal-mass 
planets. 

Whereas /3 typically represents a variable user-defined 
proxy for initial separations, the values of b and c are con¬ 
sidered to be fixed and are computed only after a suite of 
A-body simulations has finished running. Both the number 
of planets as well as the planet mass / star mass ratio under 
the assumption of initially circular orbits sets both b and c; 
these values might also be weakly dependent on ai and on 
/3. In fact, the LHS of Equation © is more often expressed 
as an unsealed time. 

The wide applicability of Equation © for systems with 
at least three planets is useful for our purposes. Although 
that equation does not let us directly relate PLATO sensi¬ 
tivities with planetary mass, as in the two-planet case, the 
relation provides an indirect pathway to do so depending on 
the architecture a reader might consider. 


3 LINK WITH PLATO 

Assume that PLATO will be able to constrain a star’s age 
to AfpLATO, such that PLATO will be able to distinguish 
(tMs/AtpLATo) time bins over a star’s main sequence life¬ 
time fins- At the level of approximation we are aiming for, 
we need not require that the number of time bins is an inte¬ 
ger, and may also treat AtpLATO as constant in time. With 
these assumptions, we can then write the number of initial 
inner planetary orbits whose total duration encompasses one 
time bin as 


3?min 


%/ GM* AtpLATO 
27T0^ 2 


( 3 ) 


One can consider x m in to be the minimum number of orbits 
needed to resolve a trend of planet frequency with stellar 
age. 

However, no resolution is possible if the planets (i) fail 
to become unstable, or (ii) all become unstable too early. In 
the first case, instability must occur before the start of the 
last time bin, such that we require 


Xuns < ( - l) X m in- (4) 

\ AC PLATO J 

This condition alone may be sufficient to detect a trend of 
decreasing planet mass with stellar age. If, however, the in¬ 
stability occurs almost exclusively within the first time bin, 
then equation © is not restrictive enough. In this case, we 
should also impose the condition 


Xuns 


> X 


min • 


( 5 ) 


We choose to leave AfpLATO as an independent 
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variable to allow for the possibility of variations in 
the stellar constraints on a star-by-star basis. Further, 
we consider 9 different main sequence stellar masses 
(0.7, 0.8, 0.9, 1.0, 1.1,1.2, 1.3,1.4, 1.5Mq), corresponding to 
the range in which asteroseismology with PLATO can be 
effective at constraining stellar ages. Because we are inter¬ 
ested only in Lagrange instability on the main sequence 
0, we must compute the main sequence lif etime 1 ms of 
these stars. To do so, we use the SSE code JHurlev et al.1 
l200dh and assume Solar metallicity. We obtain 1 ms = 
10948,7646,5607,4237,3355, and 2716 Myr for M* = 
1.0,1.1,1.2,1.3,1.4, 1.5Mq, respectively. The lower mass 
stars are on the main sequence for a duration greater than 
the current age of the Universe. Consequently, for added 
clarity, we do not draw curves in the figures of this paper 
for the 0.8Mg and 0.9 Mq tracks, as they are so similar to 
the 0.7 Mq track. 


2 planets: Minimum mass ratio for detectable trend 



3.1 Two-planet systems 

First consider equation in the context of two-planet sys¬ 
tems. Combining that with equation © yields the following 
specific condition on the planet/star mass ratio as a function 
Of AtpLATO 


Mj/Mq 


> 0.026 


1 + 0.10 log 10 


( M A\ 

( ai 

( AtpLATO ^ 

\Mq) 

V 0.2 au / 

l IGyr ) 


x 


/ £ms 

V AtpLATO 



- 5.55 


( 6 ) 

The upper panel of Figure [T] illustrates the minimum 
planet masses for which PLATO will be able to detect a 
trend of decreasing planet frequency with main sequence 
age, provided that PLATO can discover enough planets at 
these given masses to build up a large-enough sample. From 
the PLATO primary sample of stars, we expect many thou¬ 
sands of terrestrial planets with orbital periods up to about 
one year to be detectable. For a significant proportion of 
these planets, stellar activity signals will complicate the de¬ 
termination of their masses. However, by just crudely bin¬ 
ning observations, PLATO will be able to recover many of 
these systems. Disentangli ng stel l ar activity sig nals is an on¬ 
going area of research fe.g. lBastien et al.ll2014l) that will, no 
doubt, further improve the situation over the next decade. 

The first time bin is particularly important, and will 


2 Nearly all of the (transiting) planets assumed here will 
be eng ulfed by the sta r duri n g post-main-sequenc e evo¬ 
lution jMustill &: Vil laverl 120121 : iNordhaus &; Spiegell 120131 : 
lYillaver et al. I I2014I) . If the planets were further away from the 
star, then Lagrange instability could occur during e ither the gi¬ 
ant b ranch or white dwarf phases of stellar evolution (IVeras et al.l 

EH. 



AtpiATO /Gyr 


Figure 1. The minimum planet-star mass ratio [i for which 
PLATO can detect a decreasing trend of planet frequency ver¬ 
sus time for packed, Hill-stable two-planet systems. The ;r-axis 
refers to the (variable) magnitude of the stellar age constraints 
PLATO may provide. In the lower panel, the first time bin is ex¬ 
cluded, emphasizing its importance at detecting potential trends 
(shaded regions only). If ages are constrained to within 1 Gyr, 
then a trend should be detectable for planets at least as massive 
as 10 -2 Mj. 


contain a large sample of instabilities with which to com¬ 
pare. Consider the consequences of excluding that time bin; 
the lower panel of Figure [T] also imposes the restriction of 
Equation ©. If PLATO were to look for trends of planet fre¬ 
quency with stellar age only between the first and last time 
bins, then the planet mass range for which this technique 
would be successful is less than an order of magnitude. 

In any case, if we assume that PLATO will constrain 
stellar ages to within 10 per cent for stars of spectral type K7 
or earlier, then we can obtain some broad estimates for the 
planet masses for which a decreasing frequency would be de¬ 
tectable. Consider our model stars for which M+ ^ 1.0 Mq. 
Their values of £ms imply AtpLATO 1 Gyr always. There¬ 
fore, in the equal-mass case (upper panel), PLATO should 
be able to detect a decreasing frequency trend for planet 
masses ~ 10 -2 Mj and larger. This statement is also true 
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>3 planets: Necessary empirical relation for detectable trend 



AtpLATO / Gyr 



Figure 2. Values of the empirical planet instability relation ( b(3-\- 
c) for systems with at least three planets for which PLATO will 
detect a decreasing trend of planet frequency versus time. The 
solid and dashed lines coloured lines in the upper panel have 
the same designations as those in the upper panel of Figure [T] 
In the lower panel, the first time bin is excluded, such that a 
detectable trend is possible only in the shaded regions. The ± 
symbols in the y-axis labels indicate that the demarcated regions 
may be smeared out in an unpredictable fashion. If the packing 
parameter (3 is too great, then the systems will not produce a 
detectable trend. 


for 0.7 Mq ^ M* < 1.0M© because of the flatness of the 
lowest-mass curve. The bottom panel of Figure [2 also sug¬ 
gests that the trend is detectable for planet masses as low 
as those found in super-Earths. PLATO will be able to 
detect these bodies, as well as even sub-Earth planetary 
radii, and thereby improve our understanding of the cur¬ 
rently poorly-constrained radius and mass distributions of 
the lowest mass planet^. Determining masses for the light- 


3 One exception is the current record holder for the lowest mass 
planet, PSR B1257+12A (sometimes known as PSR 1257+12b), 
with a well-constrained mass of about O.O2M0 dWolszczanl 
I1994T), that resides in th e first confirmed exoplanetary system 
llWolszczan Sz Fraillll992h . 


est pl anets is difficult. Recent stu dies, e.g. IWu fe Lithwickl 
d2013l) and IWeiss fe Marcvl d2014f) . have shown potentially 
tantalising correlations but remain somewhat controversial 
due to the lack of directly-measured planetary masses. 


3.2 Systems with three or more planets 

For more complex systems with higher numbers of planets, 
we cannot make as direct a link with planetary mass as in 
the two-planet case. Instead, the key constraint is on the 
quantity ( b/3 + c). By using equation (J4J) we obtain 


6/3 + c < 10.05 + log 10 


( M *y i 

f ai r 1 1 

( AtpLATO 


} 0.2 au J 

l IGyr ) 


(7) 

We emphasize that the applicability of Equation 0 is 
limited. The intrinsic scatter in system stability causes so- 
far-unpr edictable deviation s from this formula. For example, 
Fig. 2 of lSmith fc Lissaueil (l2009l l indicates that for a given 
set of (6, c), unstable systems can occur for different values 
of /3 which vary by more than unity. 

We plot the curves resulting from Equation 0 in the 
upper panel of Fig. [5] The //-axes labels deliberately contain 
the symbol ± in order to indicate the uncertainty described 
in the last paragraph. Equation 0 implies that if the system 
is too widely separated initially (with p too large), then the 
instability will never occur quickly enough for a trend to be 
detectable. Alternatively, P can never be too small, because 
the resulting quick instabilities will be included in the first 
time bin. 

If we exclude this time bin, then the result is the bot¬ 
tom panel of Fig. [5] In this panel, the quantity (6/3 + c) is 
restricted to just about a few; too small of a range to even 
encompass the different fits from multiple authors to insta¬ 
bility ti mescales in a sim ilar region of parameter space (see 
Fig. 3 of iPu fe Wull2015ll . 

Nevertheless, and despite large phase space of systems 
with three or more planets, we can make some concret e 
statements. Consider the study of lSmith fe Lissaueil (l2009lf . 
whose equation 6 is in the same form as our Equation m and 
who derive sets of (6, c) based on extensive suites of TV-body 
numerical integrations. They found that for 5-planet sys¬ 
tems where y = 3.0035 x 10~ 6 (corresponding to Earth-mass 
planets orbiting a Solar-mass star), (6, c) = (1.012, —1.686) 
holds for 3.4 < P < 8.4. These values correspond to 
(6/3 + c) ~ (1.8, 6.8}, lower than the ranges presented in 
Fig. [5] Consequently, instabilities produced from that setup 
wou ld occur within th e first t ime bin. 

ISmith fc Lissaueil (j2009l l also considered many other 
system types; for systems with three Earth-mass planets, 
they fit values of (6, c) = (1.496, —3.142) for /3 is 3.0, which 
would give (bp + c ) « {10,11,12} for P ss {8.78,9.45, 10.12}. 
These values of /3 would place the systems in other time 
bins (besides the first). Comparing their 5-planet and 3- 
planet fits in their fig. 1 reveals that the functional form 
of the fit might change for /3 > 8.5. Quantifying how, even 
for just the equal-mass case, is difficult because of compu¬ 
tational limitations. Using IV-body simulations to predict 
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the very long term (~Gyrs) evolution of transiting planets 
(such as those which will be discovered by PLATO) pose 
an additional challenge because their semimajor axes would 
be on the scale of hund redths or tenths of an au (whereas 
ISmith fe Lissauerll2009l adopted ai = 1 au). 


4 SUMMARY 

Compared to previous experiments, the PLATO space mis¬ 
sion will significantly improve robust stellar age estimates 
by developing separate internal mass distributions for thou¬ 
sands of individual stars. Consequently, PLATO will un¬ 
veil currently undetectable trends in planetary systems. In 
this work, we demonstrated one such potential trend: de¬ 
creasing planetary frequency with main sequence age for 
marginally unstable equal-mass multi-planet systems. For 
two-planet systems, this trend could be most secure for plan¬ 
etary masses which are on the order of ten Earth masses, well 
within PLATO's planet detection capabilities. Detection of 
this trend could als o help constrain formation mechanism s 
and test the theory (lPu fe Wull2015l : IVolk fe Glad^3l2015ll 
that tightly packed inner planets were more prevalent in the 
initial stages of the lifetime of planetary systems. 
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